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Abstract 

We construct harmonic functions on random graphs given by Delaunay triangulations of ergodic 
point processes as the limit of the zero-temperature harness process. 
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1. Introduction 

Let S be an ergodic point process on M. d with intensity 1 and 5° its Palm version. Call T and 
£ the probability and expectation associated to S and 5° (we think that S and 5° are defined on 
a common probability space). The Voronoi cell of a point s in 5° is the set of sites in W 1 that 
are closer to s than to any other point in 5°. Two points are neighbors if the intersection of the 
closure of the respective Voronoi cells has dimension d — 1 . The graph with vertices 5° and edges 
given by pairs of neighbors is called the Delaunay triangulation of 5°. The goal is to construct 
a function H : 5° — >■ W 1 such that the graph with vertices H(S°) and edges {(H(s),H(s')), s 
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and s' are neighbors} has the following properties: (a) each vertex H(s) is in the barycenter of 
its neighbors and (b) \H(s) —s\/\s\ vanishes as |j| grows to infinity along any straight line. If 
such an H exists, the resulting graph is the harmonic deformation of the Delaunay triangulation 
of 5°. The search of such H has been proposed by Biskup and Berger [5], who proved its 
existence in the graph induced by the supercritical percolation cluster in Z ; their approach was 
the motivation of this paper. The harmonic function H was tacitly present in Sidoravicius and 



Sznitman [22] and in Matthieu and Piatnitski 12 ill : the function H(s) —5 is called corrector. See 
also Caputo, Faggionato and Prescott [9] for a percolation-type graph in point processes on W l . 




Figure 1 : Delaunay triangulation of a Poisson process and its harmonic deformation. The star indicates the origin 
(left) and the point H(0) (right). 



The functions from 5° to K are called surfaces. The coordinates h\, . . . , of H are harmonic 
surfaces; that is hj(s) is the average of {hj(s'), s' neighbor of s}. The sublinearity of the corrector, 
requirement (b) above, amounts to ask that hi have tilt eu the i-th canonical vector of W l . Roughly 
speaking, a surface / has tilt u (a unit vector) if (f(Ku) — Ku -u)/K converges to zero as K goes 



to iinfinity for every u e M. d (see II61I8LI15I1) 



Fixing a direction u, we construct a harmonic surface h with tilt u as the limit (and a fixed 



point) of a stochastic process introduced by Hammersley called the harness process, IU41 I18I1. The 
process is easily described by associating to each point s of S° a one-dimensional homogeneous 
Poisson process of rate 1. Fix an initial surface T|o and for each point s at the epochs x of 
the Poisson process associated to s update v\ z (s) to the average of the heights {r\i-(s'), s' is a 
neighbor of s}. It is clear that if h is harmonic, then h is invariant for this dynamics. We start 
the harness process with r|o = Y, the hyperplane defined by y(s) = Si, the i-th coordinate of s and 
show that %(•) — T|f(0) converges to h in /^(^ x P), where T is the law of the point configuration 
S° and P is the law of the dynamics. 

We prove that the tilt is invariant for the harness process for each t and in the limit when 
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t — > oo. In a finite graph the average of the square of the height differences of neighbors is 
decreasing with time for the harness process. Since essentially the same happens in infinite 
volume, the gradients of the surface converge under the harness dynamics. It remains to show 
that: (1) the limit of the gradients is a gradient field and (2) the limit is harmonic. Both statements 
follow from almost sure convergence along subsequences. 

A key ingredient of the approach is the expression of the tilt of a surface as the scalar product 
of the gradient of the surface with a specific field (see Section H]). This implies that the limiting 
surface has the same tilt as the initial one. 

2. Preliminaries and main result 

Point processes and harmonic surfaces Let S be an ergodic point process on W 1 with inten- 
sity 1; call <P its law and T, the associated expectation. The process S takes values in fA£ , the space 
of locally finite point configurations of M. d ; we use the notation s for point configurations in 9\£ 
and S for random point processes in 9{_ . The elements s of s are called points and the elements x of 
R. d are called sites. In the same way we use 3V£ ° for the space of configurations in fA£ with a point 
at the origin and s° for point configurations in that space. Let S° denote the Palm version of S. We 
can think of S° as S conditioned to have a point in the origin. If S is Poisson, then 5° = S U {0}. 
We abuse the notation and use T and £ to denote the law of 5° and its associated expectation. 
For s E fA£ let the Voronoi cell of s E s be defined by Vor(s) = {x E IR d : \x — s\ < \x — s' | , for all 
s' Es \ {s}}. If the intersection of the Voronoi cells of s and s' is a (d — 1) -dimensional surface, 
we say that s and s' are Voronoi neighbors. We consider the random graph with vertices s and 
edges {(.v,.?') : s and s' are Voronoi neighbors in s}. If S° is the Palm version of a Poisson pro- 
cess, the graph is a triangulation a.s. called the Delaunay triangulation of S°. To a site xeM^we 
associate the center of the Voronoi cell containing x: Cen(jc) = Cen(x,s) = s E s if x E Vorfs); if 
x belongs to the Voronoi cell of more than one point, use lexicographic order of the coordinates 
(or any other rule) to decide who is the center. Let 

Si := {(s,s)EM. d :sEs} 

Z 2 : = {(j,/,s)eR''xl''x^:j,/6s}. 

Functions r] : Si — > R are called surfaces and functions ^ : E2 — > M. are called fields. Denote by 
x x the translation operator: for x in R d , x x s := {s — x : s E s}. If T] (5, s) = T] (0, % s s) for every s E s 
we say that r\ is a translation invariant surface. A field ^ is covariant if C,(s' — s,s" — s,x s s) — 
C,(s',s",s) for all s, s 1 , s" E s. A field t, is a flux if C,(s,s',s) = —C,(s',s,s) for all s,s' E s. The 
conductances induced by s is the field a defined by 

a(s,s',s) := l{s and s' are Voronoi neighbors in s}. (2.1) 

The Laplacian operator is defined on surfaces r\ by 

An_(.y,s) = s)[ti(/,s)- ti(j,s)] (2.2) 

.v'es 
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The gradient of a surface T] is the field defined by 



Vti 0, s', s) = a(s, s', s) [ri (V, s) - T| (5, s)] . 
For fields C, : E2 — > M the divergence divC, : E\ — ?- M. is given by 

divC(*,s) = £a(s/,s)C(s/,s). 

Hence At] = divVr]. To simplify notation we may drop the dependence on the point configura- 
tion when it is clear from the context. The Laplacian, gradient and divergence depend on the 
conductances, but we drop this dependence in the notation, as they are fixed by (12.11) along the 
paper. 

A surface h is called harmonic for s £ 9\£ if Ah(s, s) = for all s G s. 

Pointwise tilt We say that for s G 91 a surface r\ has tilt i(r\,s) = (i ei (r\,s), . . . , i ed (r\,s)) if 
for each u G {e\, . . . the following limits for K — > ±00 exist, coincide and do not depend 

on x e R d 

i u { % *):= lim ^(Cen(x + ^),s)-T 1 (Cen W ,s)^ ^ 

Harness process Given a surface T|, let M s r\ be the surface obtained by substituting the height 
r\(s) with the average of the heights at the neighbors of s: 

{— Y a(s,s')r\(s') if s' = 5, 
<*)k (2-4) 
Tl(/) if S'^S, 

where ^(5) = E.s'es^l*^ s')- Take a point configuration s and define the generator 

L s /(ti) = ^[/(M,ti)-/(ti)]. (2.5) 

That is, at rate 1, the surface height at s is updated to the average of the heights at the neighbors 
of s. We construct this process as a function of a family of independent one-dimensional Poisson 
processes T = (T n , n = 1,2. . .) with law P. Take an arbitrary enumeration of the points, s = 
(s n ,n > 1) (for instance, s n may be the n-th closest point to the origin) and update the surface at 
s n at the epochs of T n . When the point configuration is random, say 5°, ask T to be independent 
of 5° and define the process as above to obtain a process (r\tj > 0) as a function of (5°, T), with 
the product law <£ x P, and T|o. The resulting noiseless harness process is Markov on the space 
of surfaces with generator Ls°. See Section |5] for a rigurous construction. 

Assumptions We assume that S is a stationary point process in M. d with Palm version 5°, 
satisfying the following: 

Al The law of 5 is mixing. 
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A2 For every ballfi cR d , <P(\SndB\ < d + 2) = 1. 

A3 £ exp((3a(0,5°)) < °° for some positive constant p\ The number of neighbors of the origin 
has a finite positive exponential moment. 

A4 £[(4_i(aVor(0,5°))) 2 ] <oo. The d—\ Lebesgue measure of the boundary of the Voronoi 
cell of the origin has finite second moment. 

A5 £ [£seS° a(0,s) \s\ r ] < °° for some r > 4. 
A6 T (S is periodic) = 0. 



All these assumptions are satisfied if 5 is a homogeneous Poisson process. Assumption Al 
guarantees "one dimensional" ergodicity as in (14.131) later. Assumption A2 is sufficient to define 
the Delaunay triangulation. Notice that A4 implies that the volume of the Voronoi cell of the 
origin has finite second moment: £ [(£ rf (Vor(0,S°))) 2 ] < °°. 

Assumption A6 is used on the one hand in the Appendix to identify the motion of a ran- 
dom walk on the Delaunay triangulation with the motion of the enviroment as seen from the 
walker. On the other hand ergodicity and aperiodicity of the point process imply that there exist 
measurable functions s n : 3\£ — > ~§t d such that 

Bl S- n (T Sn S°) = S n (S°), 

B2 S° = {s n (S°);neZ}, and 

B3 T Jn (5°)S° has the same distribution as 5° for every n G Z. 



This is used to extend the properties of 5° to x s S°, for all s G 5°. The point is that x s S° has the 
same law as 5° only if s is correctly chosen as was shown in [|13L l20|l for Poisson processes and 
by Timar ||23n under the condition that S is ergodic and 2>-a.s. aperiodic; see Heveling and Last 

£lS. 



Theorem 2.1. Let S° be the Palm version of the stationary point process satisfying A1-A6 and 
let y be a surface with covariant gradient, tilt I (j) G M. d and C(\Vy\ r ) < 00 for some r > 4. Then: 
(a) There exists a harmonic surface h with /i(0,S°) = and I (h) = l(y) (P-a.s. (b) ifr\ t is the 
harness process with initial condition J, then, 



lim -%(0) -h(s„) 



0, 



(2.6) 



for any nfZ, with s n as in B1-B3. (c) In dimensions d 
surface with covariant gradient and tilt 



1 and d = 2, h is the only harmonic 



Let c G M. d ; the hyperplane y(s, S°) = c ■ s, s G S° satisfies the hypotheses of the theorem with 
/(y) = c. Items (a) and (b) of the theorem say that a surface with tilt c evolving along the harness 
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process and seen from the height at the origin converges in /^(fP x P) to a harmonic surface h 
with the same tilt and with h(0) = 0. 

Let H — (h\ , . . . , hd), where /z,- is the harmonic surface obtained in Theorem l2.1l for the tilt ei. 
The graph with vertices H(S°) = (H(s), s E S°) and conductances d(H(s),H(s')) := a(s,s') is 
harmonic: 

#C0 = 4r £ a(s,s')H(s') (2.7) 
a[s) s , eS o 

that is, each point is in the barycenter of its neighbors in the neighborhood structure induced by 
the Delaunay triangulation of 5°. This graph, called the harmonic deformation of the Delaunay 
triangulation, does not coincide with the Delaunay triangulation of H(S°). 

Random walks in random graphs and martingales. Let Y t = Y t s , be the random walk on 5° 
which jumps from s to s' at rate a(s,s'). Since H(S°) is harmonic, the random walk H(Y t ) 
on H(S°) is a martingale and it satisfies the conditions of the martingale central limit theorem 
(Durrett, 11121 page 417]). So, the invariance principle holds for H(Y t ). The extension of the 
invariance principle from the walk H(Y t ) to the walk Y t requires the sublinearity in \s\ of the 
corrector %(s) = H(s)— s. 

Corrector. Mathieu and Piatnitski rf2lll and Berger and Biskup ^\ construct the corrector for 
the graph induced by the supercritical percolation cluster in Z d . Both papers prove sharp bounds 
on the asymptotic behavior of the corrector and, as a consequence, the quenched invariance 
principle for Y t for every dimension d > 2. Key ingredients in those proofs are heat kernels 
estimates obtained by Barlow yj] (in Q30 they are used just for d > 3). Sidoravicius and Sznitman 



112211 also used the corrector to obtain the quenched invariant principle for d > 4. Several papers 
obtain generalizations of similar results on subgraphs of Z d yL Q, Ej]] . Caputo, Faggionato and 
Prescott Q90 use the corrector to prove a quenched invariance principle for random walks on 
random graphs with vertices in an ergodic point process on R. d and conductances governed by 
i.i.d. energy marks. 

Uniqueness. The uniqueness of (the gradients of) a harmonic function with a given tilt has 
been proved by Biskup and Spohn [8D for the graph with conductances associtated to the bonds 
of Z rf under "ellipticity conditions" (see (5.1) and Section 5.2 in that paper) and by Biskup and 
Prescott [7] in the bond percolation setting in 7h d using "heat kernel" estimates, see Section [7] 
later. 

We obtain harmonic surfaces as limits of the zero temperature harness process. The tilt of a 
surface is obtained as a scalar product with a specific field and it is invariant for the process. This 
allows us to show that the harmonic limits have the same tilt as the initial surface. 

The paper is organized as follows. In Section[3]we give basic definitions, define the space H 
of fields as a Hilbert space and show a useful integration by parts formula. In Section|4]we show 
that the coordinates of the tilt of a surface can be seen as the inner product of its gradient with 
a specific field in M . In Section [5] we describe the Harris graphical construction of the Harness 
process. In Section [6] we prove the main theorem. Section [7] deals with the uniqueness of the 
harmonic surface in d = 2. 
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3. Point processes, fields and gradients 



Let 3\£ = 9{, (W 1 ) be the set of all locally finite subsets of R d , that is, for all s e $1, |s P\B\, 
the number of points in s P\B, is finite for every bounded set B C M. d . We consider the a-algebra 
!B(fA£), the smallest a-algebra containing the sets {s G i\£ : \sf)B\ — k}, where B is a bounded 
Borel set of M. d and k is a positive integer. 

Cesaro means and the space . Let C be the measure in £2 defined on ^ : S2 — ► R by 

C(0 = / C<*C = ^ E fl(0,5,5°)C(0,5,S°) (3.1) 

This measure is absolutely continuous with respect to the second order Campbell measure asso- 
ciated to V with density Z(w,v,s) = a{u, v,s)So(w). The space H := L2(S2,M, C) is Hilbert with 
inner product cit,- C,'), where the field (C, ■ C,') is defined by 

&-Z)(s, S ',S°)=a(s, S ',S°K( S ,s',S°)Z( S ,s f ,S°). 

If two fields C, and C,' coincide in the pairs (0, s) for all s neighbor of the origin, then their differ- 
ence has zero C -measure and hence a field in 9-( is characterized by its values at ((0,5), s neighbor 
of the origin). Define the equivalence relation C, ~ C,' if and only if ^(0,s,s°) = ^'(0,5,8°), for all 
neighbor s of the origin. Each class of equivalence in has a canonical covariant representant 
obtained by C,(s, s', s) := ^(0, s' — s, % s s) for 5, s' G s. So hereafter, when we refer to a field in H , 
we assume that it is the covariant representant. 



The space !H was previously considered by Mathieu and Piatnitski 112 lh when (5°, a) are given 
by the infinite cluster for supercritical percolation in II 1 . The Hilbert structure of this space is 
useful to obtain weak convergence for the dynamics. 

Define the Cesaro limit of a field C, : E2 — > M. by 

C(0:=lim-1- £ a( S ,s',S)^s,s',S), (3.2) 
a/>r«2|A| {v ,p A _, 

where A = A(K) := [— K,K] d C M. d . Since S is ergodic, the Pointwise Ergodic Theorem ifiol pp. 
318] implies C(Q = C(Q, fP-a.s. Analogously, for translation invariant surfaces T| we define its 
Cesaro mean C(r\) (with a slight abuse of notation) and we have C(r|) = c(r\) = £ (r|(0,S o )). 

Lemma 3.1 (Mass Transport Principle SII, 17, 20]). Let C, : E2 — > E be a covariant field such 



that either C, is nonnegative or £ Y*seS° s, S°) \ < °°. Then 

£ £C(0,*,S°) = !E £C(j,0,S°). (3.3) 

seS° seS° 

Proof. Let s„ be the maps introduced in B1-B3. Use B2 and Fubini in the first identity and 
covariance of C, in the second one to obtain 

2; £C(o,*,s°) = £^((U 3 (s°),s°) = £^(-s n (s°)A\ l{ s»)S°) 

= £2;C(^„(x,„ (S o ) s o ),0,x, i(5 o ) s o ) = £ £ C(^(5°),o,5°) = £ £C(*,o,n 
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where we used Bl in the third identity, B3 in the fourth one and Fubini and B2 again in the fifth 
one. □ 

Lemma 3.2 (Integration by parts formula). Let C,E H be a flux and § be a translation invariant 
surface satisfying £ [a(0)(|) (0)] < °°. Then 

C(V<K) = -C(<MivQ. (3.4) 

Proof Note that 

C(V<K) = U £ a(0,5,5°)V^(0,5,S o )C(0,5,5°) 

= U Y,<0,s,S°ms,S°K(0,s,S o )-U £ fl(0,*,S°)*(0,S o )C(0,j,S o ) 
= ;U £ a(0,5,5°)^(.,S o )C(0,5,5°) - ^£ [(^(0,5 )dK(0,5 )]. 

Since £ and a are covariant and § is translation invariant, a(s, s\ S )^)(s' 1 S°)C,(s, s\ S°) is covariant 
and Lemma I3T1 implies 

£ £a(0,5,5°H(5,5 o )C(0,5,5°) = £ £a(5,0,S>(0,5°)^,0,5°) 
= -£ £a(0,s,S>(0,S°K(0,s,S°) = -£[(K0,S°)divC(0,S°)]. 

We used that £ is a flux and a is symmetric. □ 



4. Tilt 

We define here the "integrated tilt" J7 (r|) for surfaces T| with covariant gradient Vr] G i# . The 
coordinates of 3 (t|) are defined as the inner product of the gradient field Vr] with a conveniently 
chosen field. We then prove that the pointwise tilt / (r|, 5°) coincides with J7 (t|), £-a.s. 

Take a unit vector u and a point configuration s°. For neighbors 5 of the origin, let b(0,s,s°) 
be the (d — 1) -dimensional side in common of the Voronoi cells of and s and let b u (Q,s,s°) be 
the projection of Z?(0, s, s°) over the hyperplane perpendicular to w, see Figure[2] Define the field 
co M by 

Q)„(0,5,s°) :=sg(5- M )a(0,5,s )4-i(^(0,^s )). (4.1) 

where i^-l 1S the (J — 1) -dimensional Lebesgue measure. By assumption A4, oo„ G 0-C and since 
Vr] is also in 0-C , we can define 

J«0l) :=C(VTi-ay) and j?(ri) := (r|), . . . , j^(r|)). (4.2) 
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Figure 2: Definition of the field co„ for u = e\ . 



Proposition 4.1. Let T] be a surface with covariant Vt) G H. Then 

i(t\,S°) = j(t\), <£ -almost surely. (4.3) 

Before proving the proposition we show a technical lemma. Let O u be the d — 1 dimensional 
hyperplane orthogonal to u: O u = {y 6 R : y • w = 0}. 

For v G O u let Z„(y) = {y + a«;a G K}, the line containing y with direction w. Fix s G 5\£ , 
define L„(y, s) := {5 G s : Vor(s) PI l u (y) ^ 0}, the set of centers of the Voronoi cells intersecting 
l u (y). Define w : R d x S 2 {0, 1} by 

, Jl ifb(s,s',s)nl u (y)^Q; 
w(y;s,s,s) = < . 

I otherwise 

the indicator that s and s' are neighbors and its boundary intersects the line /„ (y) . Define also 
0:M rf xSi ^Rby 

9(y;s,s) = ^/a + (^,5 / ,s)w(y;5,/,s), 

where <3 + (5,5 / ,5) = a(5,5 / ,S)l{(5 / -i<) > (s-u)}. In words, for s G L u (y,s), 0(y;5,s) is the neigh- 
bor of s in the direction u such that their boundary intersects / M (y). 

For x G M rf , let x* G ;( be the projection of x over the hyperplane O u . Observe that w satisfies 

w(y; s, s', s) = w(y - x* ; s - x, s' - x, x x s) , (4.4) 

and 

Q(y;s,s) -x = Q(y-x*;s-x,x x s), (4.5) 

for all xG 

9 



Figure 3: Points of L(0,s) (red and big). The horizontal line is l y . 




(4.6) 



for all A G s(R) with 1 -dimensional Lebesgue measure t\{A) < °o. 

The random set {(s- u) : s G L M (y,S)} is the one-dimensional stationary point process ob- 
tained by projecting the points of L u (y,S) to l u (y). One can think that each point s has a weight 
C,(s, 0(y; s, S) , S) . The expression on the left of (14.61) is the average of these weights for the points 
projected over A. The expression on the right of (14.61) says that this average contributes to the 
expression as much as the Lebesgue measure of the projection over O u of the boundary between 
s and its neighbor in L to its right. 

Proof. By translation invariance we can take y = and, for simplicity we take u = e\, the other 
directions are treated analogously. In this case O u = {x G M d : xi = 0}, s - u = s\, the first 
coordinate of s and x* = (0,X2, . . . ,x^). Define 



For s G 5° such that a+ (0, s, S°) = 1 , 

{s = Q(x*;0,S°),0eL u (x\S°)} = {l u (x*) nb(0,s,S°) ^ 0} = {x* G b u (0,s,S°)}. 
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= £i(A)£ £ |C(0,J,S°)| f he { x* ; o, S o )=s} l Lu{x *, S o ) (0)dx 2 ...dx d (4.7) 

r rCO J Mr 



Hence, the integral in (14771) gives a+(0, s, S°)£ d -i(b u (0, s, S )) and 

£ £g(5,5) =£i(A)£ £a + (0,.,5°)K(0,5,5°)|4 / _ 1 (^(0,.,S o ))<oo, 

because by Assumption A4 the field £d-i(b u (0,s,S°)) is in H . With the same computation, 

£ £g(s,S) = h{A)<E £ a + (0,.,S°)C(0,5,5 o )^i(^(0,5,S°)) 

*<=S ,se5° 

= h(A)<E £ a+(0,s,S°K(0,s,S> M (0,s,S°). (4.8) 

Since a+ (s, 0, 5°) C(j, 0, 5°)© M (s, 0, 5°) = a" (0, 5, S°)C(0, 5, 5°)© M (0, 5, 5°) , Lemma O implies 

£ £g(s,S) =^i(A)i £ £ fl(0,j,S°)C(0,j,S°)(0 lt (0,j,S o ) = ^i(A)c(C-(D M ). □ 

i-es 1 ses° 

Proof of Proposition ^. 1\ Again without loosing generality we take u = e\ and u- s = s\. Since 
Vr) is covariant and Cen(-) is translation invariant (in the sense that Cen(x — z,x z s) = Cen(jc,s)), 

T\(Cen(x+Ku,S),S)-T\(Cen(x,S),S) = T\(CQn(Ku,x x S) ,x x S) -r|(Cen(0,x x S),x x S). (4.9) 

Since S is stationary, if the limit in (12.31) exists, it is independent of x G W l . Let K > 0, A^ : = 
[0,K] xW 1 - 1 and define 

5k = argmax{(s-ei) : s E L(0,S) flA^}, 5^ = 0(0;^, 5), 

that is, s% is the first point of L(0, 5) to the right of A K . Let Z# = r\ (Cen(Ku, S) , S) — T| 5) and 
observe that (see Figure |3]) 

ri(Cen(^),5)-ri(Cen(0),5)= £ (ti(9(0;j,S),S) -T](j,S))+Z z (4.10) 

.seL(o,5)nA^ 

Using (14.101) . the limit as K ->■ +<*> in (1231) reads 

J M (T1,S) = Mm I £ Vti(j,9(0;j,S),5))+ Mm ^ 

= J im ^EVTi(5,e(0;^),5))l i(0)S) Wl M (^i)+ Mm ^. (4.11) 

Since (Zk)k>o is a stationary sequence, the same arguments as in the proof of Lemma l4.2l show 
that £|Z | < c(div|VTl|£(Vor(0))) < oo, and so Z K /K ->■ almost surely as K ->■ oo. Then it 
suffices to show that the first limit in (14.1 II) converges to C (Vr) ■ C0 M ) . The sum in the first term in 
(14.1 II) can be telescoped as follows: 

K-l K-l 

£ £Vii( S ,e(0;5,5),5)l L(0;5) (5)l [w+1] (5 1 ) = £(Kx fa( S) (4.12) 

k=0 stS k=0 
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where §(S) : = Y,ses ^1 ^(0; s, S) , S) Il(o.s) ( s ) 1 [o, l] (^l ) • Since the law of S is mixing, by Birkoff 's 
ergodic theorem, for integer K, 



1" 



lim - £ ty(x ku S) = <Ei$(S)} = C(Vt\-g>u) *-a.s., (4.13) 
by Lemma \4~2\ If A" e R, the result follows from the above and 

lim -\T\(Cen(Ku,S))-T\(Cen([K]u),S)\ = 2>-a.s. 

A"— >oo A 

The same arguments work for A' < 0. □ 

5. The harness process 

Given a configuration of points s G 3\£ we construct the process T|J(-,s) : s — >• R, with initial 
condition given by a surface y(-, s) : s — > R and generator given by (12.51) . 

Graphical Construction Let T = (T n ,n= 1,2,...) be a family of independent Poisson pro- 
cesses of intensity 1, T n C R. For fixed S65\[ and an arbitrary enumeration of the points of 
s = {s\,S2i ■ ■ ■ ) we use me epochs of T n to update the heights at s n as follows. Fix t > and 
define a family (M u j ; m < 5 G s) of backward simple random walks on s starting at s G s at time 
? and jumping at the epochs in T as follows. Start with B| ^ = 5; then, if x G T n and at time x+ 
the walk is at s n (that is, Bi - = s„) then the walk chooses uniformly s', one of the neighbors of 
s n with probability ^ ^ and jumps over it, setting fif , = 5'. Those jumps are performed with 

the aid of independent uniform in [0, 1] random variables U = (£/„, k, n > 1); the variable U„ is 
used to perform the k-th jump from s n . Now consider a random set of points S° with law ?P and 
assume U, T and 5° independent. Call P and £ the probability and expectation induced by (T, 
U), let P = <P x P and call E the expectation with respect to P. Denote 

p t ( S , S ',S°,T) := P(^ j0] =/|5°,r), 

the probability that , = s' conditioned on the sigma field generated by (S°,T). Define 
r|J (5, 5°, r) as the expectation of jiBt Q j) conditioned on the sigma-field generated by (5°, T): 

r\](s,S°,T) := £ ft (jy,S°,r)7(/). (5.1) 

The T] process has initial configuration r^^S , T) = 7(5) and evolves as follows. If s = s n and 
x G T n is an epoch of T,,, then 

.v"e5° 61 / 
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and r|x(s,S , T) is updated by 

i\l( S ,S°,T) = £ £ ^f^-p,-( S '', S 'X,T)y( S ') 

s'eS°s"eS° > 

La(s 7 s ,S ) / // „ T \ /c o\ 

^"ipT^' - s ' r) <5 ' 2) 

while r^^S ,? 1 ) remains unchanged for 5 V s. That is, r^(-,S°,r) = M iS T|^_(-,5 , T). 

Lemma 5.1. Given y : Si — > R wzY/2 Vy e i?/, f/ze process r\f(-,S°, •), is we// defined T-a.s. and 
has generator given by (12.51) . 

Proof. To prove that the process is well defined we need to show that the sum on the right hand 
side of (15.11) is finite fp-a.s. Proposition |9.2| in the appendix shows that 

E|Ttf(0,s°,r)|<E £ Pt (0,s f ,S°,T)\y(s',S°)\ < tc(\vy\). 

This shows that the process is almost surely well defined at the origin. Using assumption A6b 
the result is extended to all se5°. 

The fact that r\J(-,S°, •) has generator given by (12.51) . follows from (15.21) since 5° is locally 
finite fp-a.s. □ 

We have constructed the process as a deterministic function of 5° and T, the point configu- 
ration plus the time epochs associated to the points. That is, r(J is a random surface. Let (S°,T) = 
((s n ,Tn),n>l)andXs(S°,T) = ((s n -s,T n ),n>l),fovseS°. Since p t (s , s', (S° , T) ) = p t (0, s' - 
s, x,(S°, T)), y(0, 5°) = and Vy is covariant, 

T\l(s,(S°,T))= £ Pt ( S ,s',(S ,T))y(s',S°) 

s'€S° 

= £ Pt (0,s' -s,x s (S°,T))y(s' -s,x,S°) +y(s,S°) 

s'eS° 

= £ ^(0,/,x,(S o ,r))y(5 / ,x,5 o )+y(5,S°) = r|/(0,x,(S°, T)) + y(s,S°). 

If we call 

V/ (s,(5°,r)):=riJ(0,x,(5 o ,r)), (5.3) 

then the process at time t is the sum of the translation invariant surface \|/ r and the initial condition 
y. That is, 

% = ¥f + Y (5-4) 

In particular, it follows that Vr|J is a covariant (random) field P-a.s. 

The dependence of r\] on (S°,r) will be dropped from the notation when clear from the 
context. 
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Extension of the Hilbert space H to include the randomness coming from the process We 

consider the probabilistic space where S c ',T,U are defined as independent processes and abuse 
notation by calling C the Campbell measure on E2 associated to 5°, T, U: 

c(Vr\J) := ££(Vr|J) = E(Vitf) 

The following bound -shown in the Appendix- implies that the process is well defined as an 
element in H for all time. 

Lemma 5.2. IfC(\Vy\ r ) < °° then 

C(|V ¥ Jr) < 2 r c(\Vy\ r )m r (t)<°°, 

where m r {t) denotes the r — th moment of a Poisson random variable with mean t. 

As a consequence of (15.41) and Lemma [5^21 the tilt is invariant under the dynamics: 
Proposition 5.3. For all unitary u E M. d and covariant surface y with Vy e H , 

Ju(vJ) =3 u {i). 

for all t > 0. 

Proof. First observe that with ^-probability one we have, 

divco M (0, 1 S°)= £(D B (0,j,S°) = i £ 4_i(^(0,5,5°))-^ £ £ d -i(b u (0,s,S°)) = 0, 

(s-u)>0 (s-u)<0 

because each term in the substraction is the (d — 1) -dimensional Lebesgue measure of the pro- 
jection of the Voronoi cell of the origin over the hyperplane orthogonal to u. Then, 

JuCn7) = C (VnJ • (Ou) = C (Vy- G)„) + c (V\|/ f • 0)„) 

= c(VY-(o M )-c(Vrdiva) B ) = c(Vy-co„) = j7 l( (y). 

where we used (15.41) in the second identity, the integration-by-parts Lemma l3?2l in the third iden- 
tity as V|/J is a translation invariant surface and divco ; , = in the fourth identity. □ 

6. The process converges to a harmonic surface 

In this section we show that if y is a surface with tilt l(y), whose gradient is in 9{ and has 
more than 4 moments, then there exists a surface h with V/z e H such that VtjJ converges strongly 
in ^ to V/z. Furthermore /z is harmonic and has the same tilt as y. We split the proof into several 
lemmas. 
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Lemma 6.1. IfC(\Vy\ r ) < °°for some r > A, then for all t > 

j t cml\ 2 ) = -2E[a(0)- 1 |Aii7(0,S o ,r) 



2 



(6.1) 



Proof. We drop the dependence on the initial condition y, S° and T and write V[ t = "(]](■, S°, T). 
Let <X 2 = Ui„eV2 *- ne e P ocns corresponding to sites in V 2 , the set of second neighbors of the 
origin. Define the events 

Fi := Fi(t,h) = {\<r 2 n[t,t + h]\ = l}; 
f 1jS := F 1>s (t,t+h) = F l n{\T(s)n[t,t+h]\ = l}n{ s eV 2 }; 

F 2 := F 2 (t,h) = {\T 2 n[t,t + h}\>2}. 



Given S, T 2 is a Poisson process with intensity \V 2 \, hence 



P(Fi|5°) = E[l Fl |5°] = \V 2 \he-^\ (6.2) 

¥(F 1:S \S°) = he-^-H V2 (s), (6.3) 

F(F 2 \S°) < h 2 \V 2 \ 2 . (6.4) 

We have to compute 

E£a(0,^)(|Vri m (0^)| 2 -|Vri f (0,5)| 2 )(l Fl +l F2 ) = / + // (6.5) 

seS° 

We use 

|VW0,s)| 2 -|Vr| f (0,s)| 2 = [Vri /+/7 (0,5)-Vri ? (0,5)] 2 + 2Vri f (0,5)[Vri f+/l (0,5)-Vri f (0,5)], 

to compute each term in (16.51) . Assume F\ occurs. 

• If the mark is neither at the origin nor at a neighbor of it, then a(0,s) = 0, Vr^^O,.?) = 
VT|f(0,s), and the difference is zero. 

• If the mark is at the origin and a(0, s) = 1, 

|Vri /+/7 (0,5)| 2 -|Vri ? (0^)| 2 = [-M %(0) + r| f (0)] 2 + 2Vr| / (0,s)[-Mor| / (0) +%(0)] 

= -2Vri,(0,5)A*ri f (0) + |A*ri f (0)| 2 . (6.6) 

• If the mark is at some s such that a (0,s) = 1, we have Vr^^O,.?') = Vr| f (0,s'), for all 
s' ^ 5. So 

|Vri /+/7 (0,5)| 2 -|Vri f (0,5)| 2 = [M s t\ t (s) -r\ t (s)] 2 +2VT\ t (0,s)[M s T\ t (s) 

= 2Vri f (0,s)A*r| f (s) + |A*r| ? (s)| 2 . (6.7) 
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Given 5°, the process T2 D [t,t + h] is independent of r| f , so conditioning on 5° by (16.21) . (16.31) . 
(16.61) and (16.71) . we get that the first term in (16.51) equals 



fcEfe-l^ £ a(0,5)(2Vri f (0,5)VA*ri f (0,5) + |A*ri f (5)| 2 +|A*ri ? (0)| : 

By monotone convergence, 
1 



S^ E ( E «(0,5)(|V %+/l (0,5)| 2 -|VTi r (0,5)| 2 )l Fl 



E(£a(0^)Vri f (0^)VAX(0,5))+^E(£a(0,5)(|AXW| 2 +|A*ri f (0)| 2 ) 
e( £ a(0,s)VTi f (0, j)VA*Ti f (0, j)) +E(a(0)|A*r| 



s€S° 



(6.8) 



by the Mass Transport Principle (13.31) . Let l/p + l/q= 1 and £(0,5) := \p 2 , then for any time t', 
by means of (15.41) and Lemma [5^21 the second term in (16.51) reads 



E(£|VTv(0,5)| 2 l ft 

\es° 



h 



c(\Vr\ t >\ Z Q 



h 



E(E |V^(0,5)| 2 ^)l 1/P [Ea(0)l F2 



1/? 



< (Am 2p (0+5) £|y 2 | : 



l/q 



for constants A, 5 > 0, where m r (t) is the r-th moment of a Poisson random variable with mean t. 
Choosing q < 2 and applying this bound for t' = t and t' = t + h we get 

lim^E(£a(0,5)(|VTi m (0^)| 2 -|VTi,(0^)| 2 )l F2 ) =0. (6.9) 



From (16.51 ). (16.81 ). ( 16.91) and the integration by parts formula we obtain 



d 
dt 



C(\Vr\ t \ 2 ) = 2c(Vr|,VA*r| f )+E[a(0)|A*r| / | 2 ] = -E[a(0)|A*%l 



□ 



Corollary 6.2. If y satisfies the hypotheses in Lemma \6J] then 

(a) <T ( I Vr|7| 2 ) 15 non-increasing in t; 

(b) C(\ Vr|J| 2 ) is strictly decreasing at time t if and only ifrtf is not harmonic for (a,S°); 

(c) \im t ^ oo a(0)~ 1 Ar{f(0) = 0, F-a.s. and in L2(P), 2-a.s.; 

(d) lim^ooAr)/ = F-a.s. 
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Proof. Let 

^:=*-(0)|AV(0)| 2 
a(0) 

Lemma ISTTl implies f^"E[Z t ]dt < °°. Fix to = and denote = to < t\ < ti < . . . the ordered 
epochs of the superposition of the Poisson processes associated to the point at the origin and its 
neighbors. This is a Poisson process with intensity a(0) + 1. For each n > 0, given 5°, Z tn is 
independent of (t n+ \ —t n ). Hence, 

/ EZ t dt = E Z t dt=YE[Z tk (t k+1 -t k )} = £E( * ) <co. 

Hence, 

oo 

V Aril (0) < °° and lim Ar^O) = P-a.s. 

k=0 

The L2 (P) convergence follows by dominate convergence using that Ar| J(0) < Y,k=o AiU (0) . □ 

Proof of (a) and (b) of Theorem 12.11 In the notation we drop the dependence on the initial sur- 
face y. We want to prove the existence of a harmonic surface h, with covariant V/z and such that 
for all n E Z 

lim Ea(0,5 ;j )[Vti f (0,5,0 -V/*(0,s„)] 2 = 0. 

where (s n ,n E Z) is the enumeration of 5° given in B1-B3. 
Observe that 

E\a(0,s n )(Vr\ t (0,s n )-Vh(0,s n ))\ 2 < E£ a(0,s)\VT\,(0,s) - Vh(0,s)\ 2 = 2c(\Vr\ t - Vh\ 2 ). 

seS° 

So, it is enough to show that Vr], — > Vh strongly in H . 

Existence of the limit. Since by Corollary |6.2l c(\ Vr[ f | 2 ) is bounded, V% is weakly compact, 
and hence for every sequence {tk}k>o, there exists a subsequence {tkj}j>o and a field ^ M €^ 
such that 

limc(Vm ■C) = C(Coo-0, for all Ce^"- (6.10) 

Uniqueness of the limit. Let {^}yt>o be a subsequence such that Vl|^ — ^00. 
By (1541). 

C(|^| 2 ) = limC(VTV^) = c(VyCoo)+limC(V\|/^-Coo). (6.11) 

where \j/ f is defined in (15.31) . Integrating by parts and using Holder, 

I C ( V\|/ ft ■ Coo) I = lim I C (Vy* • VT] fj ) I = lim I C (\|% • Aty ) | 

< limE(a(0)|v ) / r J 2 ) 1 / 2 E(a(0)- 1 |Ari f .| 2 ) 1 / 2 = 0, (6.12) 

;-J.co 
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by Corollary 16 .21 Therefore, 

C(|W 2 ) = c(Vy-U)- (6-13) 

Let Vt|^ — ^ and Vt|^ — ^ subsequences converging to two weak limits and C^. By (16.101 ) 
and (16.111) . 

C(Coc-CL) = lim^VTvQ = C(VY.Q + limC(V ¥ ,,-a = C(|C| 2 ), (6.14) 

by (EH and dBTDT) . The same holds for ^ and so C(\C,l\ 2 ) = C(\^\ 2 ) = C(U ■ This 
implies C(|^oo — CJ 2 ) = 0, i.e. there is a unique limit point. 

Strong convergence. By (15.41) and integration by parts, 

c(\VT} t \ 2 ) = c(VtVtiO + c(Vy,VtiO = c(V-yV%)-c(v,ATiO. (6.15) 
From Holder's inequality, 

(c(v|/ ? A^)) 2 < E(a(0)| ¥f (0)| 2 )E(^|^). (6.16) 
Since by Lemma I6TI E is integrable, there exists a subsequence (tk)k>o such that 

jt->oo V a(p) 
From Lemma [9731 in the appendix, 

££|y(5j ' )|2 , A (0) ,2 

lim ^ — t k E ] Uk \ Jl = 0. (6.17) 

k^°° tk a(0) 

Using (16.131) . (16.151) and (16.171) . C(|Vr| ; | 2 ) ->■ c(|£«>| 2 ), and hence Vr| f converges strongly in H 

tO ^oo. 

Zero divergence. By Jensen's inequality and using a(0) > 2, we get 



limE(a(0)- 2 |Ari f -div^| 2 ) < limE(a(0)" 1 Y a(0,s)(VT\ t (0,s) -^(0,s)) 2 ) 

< lim C(\Vr\ t - U| 2 ) = 0. 

f — .co 

It follows by Corollary I6.2l that 

div^oo = P-a.s. (6.18) 

Covariance. A field £ G # is characterized by its values on the edges leaving the origin. 
Therefore, by taking the covariant canonical representant defined by £00(5, s' ,S°) := ^(O,^ — 
s, l s S°), we can consider ^ to be covariant. 
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Gradient field. To show that is a gradient field we prove that it verifies the co-cycle prop- 
erty, that is there exists C 9£, with t(S° G = 1 and such that for all s G i\£ * and every 
closed path s io ,s iv ...,s ik = s io G s with a{s tj , t ) = 1 , j = 1 , . . . , Jc we have Y!}= \ C~ ( s > s ij- u s ) = 
0. 

Let /i,/?z G Z. Since 61(^)2, i' m )VT|f(5 n ,5 m ) >■ (5^, 5,^)^00(5^5 

s m ), we have a subsequence 

that converges almost surely. Denote by 9{n, m C 9£ the set where convergence holds. Using a 
standard diagonal argument we get a subsequence (V&)fc>0 sucn mat 

a(j„, ^Vt^^Jm) a(5„,5 m )Coo(^,*m) for all n.m G Z. 

Define 9£* = (\ nm ^z^Gi,m- Since the co-cycle property holds for every t the a.s. convergence 
implies the co-cycle property for ^00. 

Tilt. The tilt is a continuous functional in H and it is constant for the dynamics by Proposi- 
tion [531 Hence the limit £00 has the same tilt as the initial surface. This completes the proof of 
(a) and (b) of the theorem. □ 



7. Uniqueness of harmonic surfaces in d = 2. 

In this section we prove uniqueness (up to an additive constant) of the harmonic surface 
with covariant gradient for d = 2. Observe that in dimension one the harmonic function with 
a given tilt can be explicitly computed and hence the uniquness follows immediately. To prove 
uniqueness for d = 2 we use the following result. 



Theorem 7.1 (Theorem 5.1 of Berger and Biskup Q3D )- For c G M 2 , let y(s) = c-s, and h be a 
harmonic surface for a(-,-,S°) with covariant gradient in H and tilt I (h) = I (y). Then 

Urn- max {\h(s) - c ■ s\} = 0, fp-a.s. (7.1) 

n s eS°n[-n,n] 2 



We omit the proof; it follows |5J], details can be found in \ j l6 \. Berger and Biskup [5] use 
this theorem to show uniqueness of the harmonic surface on the supercritical bond-percolation 
cluster in Z 2 ; we adapt their proof to our case. Theorem 2.4 of Bikup and Prescott [Q] proves 
(17.11) for bond percolation in Z d for all d > 2 under "heat kernel estimates" assumptions, see 
(2.17) and (2.18) in that paper. These estimates are to be established in our setting. 

Proof of (c) of Theorem 12.11 It is enough to show that if h is a harmonic surface with I (h) = 0, 
then V/z = or, equivalently, C(|V/z| 2 ) = 0. From the considerations after (13.21 ), if V/z G 0-C then, 
with probability 1, 
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Let S n = S fl [—n, n] 2 . Using that h is harmonic rewrite the sum at the right hand side as 

£ a(s, s')\Vh(s,s') | 2 = £ VfcC*,/) " £ h ( s ) £ «0mO VA(j,/) 

ses« ses n seS n s'eS 

s'eS s'eS 

= £ a(s,s')h(s')Vh(s,s'). 

ses„ 
s'eS 

Using harmonicity again, we obtain 

^a(s,s')\Vh(s,s')\ 2 = £ a(s,s')h(s')Vh(s,s') + £ a(s,s')h(s')Vh(s,s'). 

S£S n seS n seS„ 

s'eS s'eSn s'es\s„ 

= £ a(s' ,s)h(s)Vh(s',s) + £ a(s,s')h(s')Vh(s,s'). 

s'eSn s'es\s„ 

= — £ a(5 / ,5 , )/z(5)V/z(5 / ,5') + £ a(5,5 / )/z(5 / )V/z(5,5 / ). 

*G5„ seS„ 
s'65\5 B s'eS\S n 

= £ fl(5,/)(^(5)+/z(/))V^(5,5 / ) 

s'eS\s„ 

Then, with ^-probability 1, 

C(\Vh\ 2 ) = lim V a(5,5 / )(/z(j / )+/z(^))V/z(5,/). 

Since this limit exists a.s., we are done if we can show that the r.h.s converges to zero in proba- 
bility. Observe that 

£ a(s,s')(h(s)+h(s))Vh(s,s) < max {a(s,J)\h(s) +h(s')\} £ a(s,s)\Vh(s,s)\. 
Let A„ := {There exists s E S n and s' G S\S2« such that a (s, s') = 1 }, and observe that 



<P(A n ) < £ £ £ a(*/) < £ £ £ fl (,,y)i— g-L < -£[£a(0,,)|.| 4 ] 

seSntfesySto seS n s'eS\S 2n n n seS° 



Therefore, by Borel-Cantelli, the fact that 1(h) = and Theorem |7.1[ given £ we can take n big 
enough such that 

(2ny l max {a(s,s')\h(s) + h(s')\} < - max{\h(s)\} < 6. 

seS n ,s'eS\S n nseS 2 „ 
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It follows that 



lim (2n)~ l max {a(s,J)\h(s) +h(s')\} = 0, 2>-a.s. 



and therefore it is enough to show that there exists a sequence (Z n )„>i such that 

Z„ > -US) := - £ £ a(s,s')\Vh(s,s% 

n n seS„ seS\S„ 

almost surely and Z n converges in probability. 

Given B,B> e <B(R 2 ), let ^(S) := LsesLs'es ■a{s,s , ,g)\7h(s,s , ,S)\l B (s)lff(s'), and ob- 
serve that by the refined Campbell formula and the covariance of V/i and a 

£<t>s,s'= / ' E £ a(5,/,x_,5)|V/i(5,/,x_. s 5)|l s (5)l B '(/)J5 



£ a(s, s' + 5, %- s S) | V/i(s, s' + 5, x_,S) | l fl (s) l s / (s' + s) 



s'eS 



ds 



s'eS 



ds 



<E £ a{0,s',S)\Vh(0,s\S)\£{Bnx s ,B') 

s'es 



(7.2) 



Let S n = [— n,n] 2 and be the family of half -planes defined by the borders of B n , and 
disjoint from B n . It is clear that 

£ a{s,s')\Vh{sJ)\< £ K„ s (5). 

s'eS\S„ 

We show the convergence of ^6„,b(5) for a fixed B e X„. The convergence of the other terms 
follows from the same arguments. 

Before proceeding, we have yet another approximation to take care of. Let //„ = lx [n, +°°), 
G„ = [— n,n] x (— <»,n], and observe that 

Let us see what happens with a fixed line first. To do that, let G = [0, 1] x M~ and G° n = 
[—n,n] x R~. If we define T = x ei , by the covariance of V/z and Birkhoff Ergodic Theorem, it 
follows that 

1 1 " _1 

lim -<^ G « Bb(S) = - lim V § G , Ho (T k S) = 2-E [§g,h (S)} < °° a .s. 

ii — VI 11 1/1 n — ^nn 



k=—n 



By the covariance of V/7 it follows that 

lim2'(|^ G „ ;/f( ,(5)-2E[^ Gi//o (5)]| >en) = lim (p(\^ G o Mo (x ne2 S) -2<E [§g,h (S)}\ > en) = 0, 



71—5-°° 

and the result follows. 



□ 
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8. Final Comments 



8.1. Invariance Principle 

The key ingredient to obtain an invariance principle from the existence of a harmonic defor- 
mation of the original graph is a uniform sublinear bound of the corrector as in (17.11) . GrisiJ 1611 
obtained this bound for the Poisson process following the arguments of Berger and Biskup pD in 
d = 2. Hence the quenched invariance principle holds in the Delaunay triangulation of a Poisson 
process. Presumably this also holds for a non-periodic ergodic process satisfying assumptions 
A1-A7. For d > 3 the proofs of a quenched invariance principle in the percolation setting and 
related models rely on heat kernel estimates like those obtained by Barlow [1], which do not 
follow from the sublinear behavior of the corrector along lines. An extension of these bounds to 
our case are to be obtained. 



8.2. The process trajectory is orthogonal to the space of harmonic surfaces. 

Since the tilt in the direction u G M. d is a continuous functional in H , by Riesz Theorem, there 
exist a field 0)„ G itf such that the tilt is given by the scalar product with (D u . In our case, we have 
found explicitly that field (the one given in (14.11) ). 

Given an initial condition y, the process = T] ] — y is a translation invariant surface and has 
zero tilt. The convergence of V\j/ ? follows from the convergence of Vr|/, and the limiting field is 
the gradient of the corrector V%y := V/z — Vy, for h given by Theorem 12.11 Integrating by parts 
and using translation invariance, for £ G H with divC, = 0, 

C (V(y - TiJ) Q = - C ( (y - TlJ)divQ =0, for all t > 0. 

Hence y — r\] is orthogonal to the subspace of fields in H with zero divergence (which contains 
the gradients of all harmonic surfaces). In fact, V/z is the orthogonal projection of Vy over this 
subspace. In particular, we have 

Vy= V/z + (Vy- Vh) (8.1) 

Mathieu and Piatnitski rf2~lll also consider £2(22, C). Equation (18.11 ) corresponds to their 
decomposition of the space as £2(22, ) = ^2° Z ©^2°'- Taking y,-^) := ■ s), i < d, the surface 
X := (Xyi ? • • • 1 Xy d ) 1S wnat m ey call the corrector. 

8.3. Regularization effect. 

The regularization effect observed in Fig. [T]can be explicitly formulated as follows. If one 
takes n arbitrary points si,...,s n G M 2 , the barycenter minimizes the following sum of scalar 
products 

n y 11 

argmin V [(s k -x) ■ (s k+l -x)] = -Ys k . (8.2) 
*^ 2 k=\ n k=\ 
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where s n +\ =s\. Take a point configuration s and s,s' Es neighbors in the Delaunay triangulation 
of s. The directed edge (s,s r ) is shared by the triangles ss'a + and ss'<X-, where a + (s,s') is the 
first common neighbor of s and s' in the clockwise direction from s' — s and a_ (s, s') is the other 
common neighbor. We show the following extension of (18.21) to harmonic surfaces. 

Lemma 8.1. Let S be a stationary point process. Then the harmonic deformation of the Delaunay 
triangulation ofS minimizes 

■E £a(0,s)[G(s)-G(a+(0,s))] (8.3) 

among deformations G : 5° t— > W 1 ofS such that G(0) = and the corrector G(s) — s has coor- 
dinates with gradient in M . 

We prove this Lemma below. Given a surface T] define the fields £}, : S2 — > K by 

:= a(s,s')VT\(s,a±(s,s')). 

Any two surfaces r\,§: Si — > K satisfy 

C(Vri^) = c(Cv^). (8.4) 

Also note that 

^a(s,s%l(s,s')=Ar\(s) and £ a(s,s%l(s',s) = 0. (8.5) 

If (f> is a translation invariant surface (that is (|)(j,s) = (|)(0,x s s)) then, by the mass transport prin- 
ciple, 

2c(v^l) = <l £ fl(o,j)V4>(o,5)a(o,j) 

= £(|)(0) ^a(0,5)S(^0)-£(t)(0) £ a(0,j)Cl(0,j) 

= -2c(<|>Ati)=2<:(V<|>Vti), (8.6) 

where the first identity in the bottom line follows from (18.51) and the second one by the integration 
by parts formula. 



Lemma 8.2. 



d -C(Vr\ t Q) = ^(|Vr| f | 2 ) = — E \a(oy l \Ar\ t (0,S°) 



dt ^ 2dt 

Proof Using (|84|) and V% = Vy+V\|/ f , 

c(v%B) = c(v Y ft) + c(Clv ¥f ) + c(v ¥f S) 

= C ( Vyft) + C ( VyV ¥f ) + c (V\|/ f Vr| f ) 

= C(Vy^) + c(VyVri f ) + C(V^Vri f )-c(|Vy| 2 ) 

= C(VyC + ) + c(|V % | 2 )-C(|Vy| 2 ), 



(8.7) 
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where the second identity follows from (18.61) . This shows the first identity in (18.71 ); the second 
identity is (I6TTT) . □ 



Proof of Lemma I&71 Lemma [8T2l shows that c(V%yf) is non-increasing, and that it is strictly 
decreasing if and only if r\ t is not harmonic and hence 

C (Vg^j_) =0 if and only if g is harmonic 

Taking gi as the coordinates of G and using that G(0) = 0, we get (18.31) . □ 



8.4. Some simulations. 

The first two pictures in Figure |4] show level curves of a linear interpolation of the surface 
y— h. In the first one some level curves are drawn. From blue (minimum) to red (maximum). 
The level curve of zero is drawn in green. In the second one the sublevel set of zero is drawn in 
blue and the superlevel set is drawn in red. The black curve is the level set of zero. 

The next picture is the Voronoi tesselation of the harmonic points. The Delaunay triangula- 
tion of this points does not necessarily coincide with the harmonic deformation of the original 
Delaunay triangulation. It is easy to construct examples where this in fact happens, and it can be 
seen in simulations. However, it can be appreciated in simulations that the density of triangles in 
the harmonic graph that are not Delaunay triangles is very low, as shown in Figure [5J Finally, on 
the bottom-right of FigurelH the level curves of the harmonic surface with tilt ( 1 , 0) is shown, that 
is the limit of the dynamics with initial condition given by the hyperplane y(x, y) = x. Observe 
that the surface is pretty close to the original condition. 



9. Appendix. The Random Walk and the Environment Process 

This appendix collects some technical results used in Section [5J The environment seen from 
the particle was used by De Masi et al. Ill ill to show the annealed invariance principle for the 
random walk in the supercritical bond-percolation cluster. We adapt some of those results to our 
setting. 

Let s G 9l and s G s. Let X* be a discrete time random walk on s with law P s defined by 
Xq = s and for n > 1 , 

P s (XS = s"\X*_ l =s>)= . 

That is, the walk starts at s and if it is at s' G s, then it chooses a neighbor uniformly at random 
and jumps over it. Let E s be the expectation with respect to P s . 

To build the continuous time walk, let N = {T^k G N} be a rate 1 homogeneous Poisson 
Process in R+, independent of (X n ) n >o, and define 

X t :=X N{t) , (9.1) 



a(s',s",s) 
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Figure 5: Delaunay triangulation of harmonic points (black and thick dashed) vs. harmonic graph (yellow) 

where N(t) = jiVD (0, t] | is the number of points of N in the interval (0, t] . Let P s = P s <g) Q, where 
Q is the law of N in (R+), S (v£ (K + )))- The law of X? coincides with the law of the walk 
B? j defined in Section [5l so that the results below hold for B® Q j . 

Given the process Xf r (with initial state s G s), define the process 

This process can be thought as the environment as seen from the particle moving according to 
X„. The process s° is Markov with values in £V^° (i.e. for all n, G s°). We use P s to denote the 
law of s° in with initial state s . 

Let M be the set of aperiodic s: 

M = {s G 9i : i x s ^ s for all x G R d , x ^ 0}. (9.2) 

If s is aperiodic, then the trajectory of s„ determines univoquely the trajectory of the walk X®. 
The Poisson Process is aperiodic almost surely. 

Let S be an ergodic point process in W 1 , with Palm version 5°. Denote by Q the probability 
measure on , S given by 

for bounded measurable / : 9\£ — >• R. 
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Lemma 9.1. The process (s°)„>o is reversible and ergodic under Q,. 

Proof. To check reversibility, let /, g : 5\£ — > K. be bounded measurable functions. Define (|>(s, 5', 5°) = 
a(5,5 / ,5 )/(x x 5 )g(x J '5°) and observe that § is covariant and integrable, and therefore, by the 
Mass Transport Principle (Lemma [37Tb 

f £ s o/(s°)s(s^(Js°) = (l/E fl (0))E £ a(0,s,S°)/(S°Mx,S°) 

= (l/E[a(0)])E £ a(0,s,S°)/(T. ? S°MS°) 

To show ergodicity, let A £ CB(9l°) be an invariant set for the dynamics, that is A is such that 
Sq £ A implies Sj £ A. This implies that for any neighbor s of the origin, t s Sq £ A. Iterating the 
argument one shows that, if s° £ A then x 5 s° £ A for every s £ s°. Therefore, 

fP (A) := fP(S° £ A) = lim -L V l Ts s° e A e {0, 1} 

A^|A| s £s° ' 

and, as Q, < fp and fP < Q,, it follows that Q, (A) £ {0,1}. □ 
Proposition 9.2. Lef r > 1 and J be a surface with covariant gradient. If c := 2c(|Vyj r ) < 00 

E\j(X t )-y(X )\ r = <£E S o\y{X t )-y{X Q )\ r < <£{a{Q)E s °\y{X t ) - y(X )\ r ) < cm r (t), 
where m r {t) is the r-th moment of a Poisson random variable with mean t. 

Proof. Suppose, without loss of generality, that y(0) = 0, and observe that 

00 00 n 

<E(E s »\y(X t )\ r ) = Y,-E(E s »\y(X n )\ r l m=n ) = £ e^I^D^- 

n=l n=l H - 

On the other hand, by Holder's inequality 

lY(*n)| r = I £ (y(Xk) -Y(^-i))l r < n r - x £ \vy(x k _ u x k )\ r . 

k=\ k=\ 

Finally, as X k — Xk-1 depends only on S£ and s#_i, by the stationarity of s n under Q, (Lemma 
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[97Tb . it follows that 

<E{E S o\y{X t )\ r ) < ^n^ 1 £<E(£ S o|VY(X fe _ 1 ,Z,,S°)r) e -^ 



n=l /t=l 

< ^n^ 1 ££(fl(0)^|VY(0,X fe -l fe _ 1 ,s,_ 1 )r) e -^ 

n=l fc=l 

= £n''£(fl(0)^|VY(0,X 1 ,5 o )r)^-r 

= £(£a(o,5)|v Y (o, s ,s )|'>^). □ 

To obtain estimates for cQVT|J| r ) we study the process of the environment as seen from the 
random walker on 5°, as in 111 ill . The law of 5° is reversible and ergodic for this process, which 
allows us to make estimates on the original random walk. Let B? ( q , as in Section [5] be a random 

walk on the points of S° starting at G 5°, and denote its law by P^ a . 

Proof of Lemma 1572] From the covariance of Vy we can assume, without loss of generality, that 
y(0) = 0. By the Mass Transport Principle Lemma I37T1 and Proposition 19.21 

C(\Vy t \ r )= ^E£ fl (0,j)|V V/ (0,j)r < T- 2 E 1 £a(0,s)i\y t (0)\ r + \Vt(0)\ r ] 

1 seS° s€S° 

< 2 r - l Ea(0)\y t (0)\ r < 2^-^(0) |y(< 0l )r < 2 r C (\Vy\ r )m r (t) fP-a.s.. □ 



The following Lemma is a part of the proof of Theorem 2. 1 in DJJJ] ; we omit the proof. 
Lemma 9.3. If a surface y satisfies 



2 - 



££a(0,,)|y(,)| 2 < 
then 

h ^igg<, ana lim *W0)Ehffll 2 ) <oo . (9 . 3) 

n— >°° n t— t 
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